Discrete charging of a quantum dot strongly coupled to external leads 
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We examine a quantum dot with A'^dot levels which is strongly coupled to leads for varying number 
of channels in the leads. It is shown both analytically and numerically that for strong couplings 
between the dot and the leads, at least A^dot — N bound states (akin to subradiant states in optics) 
remain on the dot. These bound states exhibit discrete charging and, for a significant range of 
charging energies, strong Coulomb blockade behavior as function of the chemical potential. The 
physics changes for large charging energy where the same (superradiant) state is repeatedly charged. 
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It is well known that the number of electrons in a 
weakly coupled quantum dot changes discretely as func- 
tion of the chemical potential. This phenomenon is the 
basis for the application of such dots as single electron 
transistors Q . When the coupling to the external leads is 
weak, it may be treated as a perturbation and results in 
a broadening V = 7r7Vz/|T^p of the states of the uncoupled 
dot. [u is the density of states (DOS) in a lead, TV the 
number of leads, and V is the overlap matrix element be- 
tween a state in the dot and a typical state of the leads.) 
Usually, one expects the discrete features of the dot to 
be lost once F is larger than the typical level spacing A 
in the dot. This corresponds to the requirement that the 
dimensionless conductance through the dot, gdot = T/A, 
should be larger than one. 

Electron-electron interaction in the dot results in the 
Coulomb blockade phenomenon 0]. Within the "ortho- 
dox model" 0, the chemical potential change needed to 
add an additional electron to a weakly coupled dot is no 
longer A but rather A-|-e^/C, where C is the capacitance 
of the dot. For stronger coupling to the leads, suppres- 
sion of the Coulomb blockade is predicted 01, and once 
r ^ A only a weak remnant of the discreteness of the 
dot is expected 

In this picture the only relevant condition for the the 
appearance of discrete features in the dot is g^ot < 1- 
Nevertheless, one may extrapolate from the Dicke effect 
0, 0] that a quantum dot strongly coupled to a lead will 
also show sharp resonances. In an extreme strong cou- 
pling limit, we can think of the dot levels as degenerate, 
analogous to the identical resonances of the atoms in the 
Dicke effect. These degenerate levels are coupled via the 
continuum of lead states which is akin to the coupling 
of the atoms by the radiation field. Indeed, a resonance 
in a two orbital dot strongly coupled to two leads was 
seen in Ref. 0- The tunneling DOS of a non- interacting 
two-orbital dot strongly coupled to a single lead tends 
towards a delta-like peak when gjot — > oo |^. It was 
later shown that when a dot of A^dot states is strongly 
connected to a single lead, iVdot — 1 delta-like peaks in the 



tunneling DOS remain for gjot ^ oo. The case of two 
leads connected to a dot with A^dot — 2 shows interest- 
ing dependence on the overlap matrix element Vk,i (i.e., 
the i-th level overlap with fc-th lead) lOj. For identical 
matrix elements one sharp peak in the local density of 
states remains, while if one of the four matrix elements 
has a different sign, no such features are observed. 

In this letter we show that in the limit g^ot ~> oo, the 
relevant parameter in determining the number of bound 
states in the dot for generic dot-lead coupling is the num- 
ber of channels N of the leads, or equivalently, the di- 
mensionless conductance g of the leads. We demonstrate 
that for N channels coupled to a dot, TVdot — N states 
remain bound to the dot, except when the coupling ma- 
trix elements between leads and dot are independent of 
either dot level or channel index. For the latter cases, 
-^dot — 1 states remain bound to the dot. Using numerical 
density-matrix renormalization group (DMRG) as well as 
analytical arguments, we show that these bound states 
exhibit discrete charging as well as Coulomb blockade 
except for very large charging energies. Thus, in order 
to wash out all discrete features of a dot in the limit of 
strong coupling, one or more leads of total dimension- 
less conductance g > N^ot must be connected. With 
increasing coupling strength, a dot coupled to leads will 
evolve from N^ot bound states at very weak coupling, to 
no bound states at intermediate coupling ^ and finally 
to A'^dot — g bound states at strong coupling. Criteria for 
the different regimes and experimental realizations will 
be discussed. It is interesting to note that the physics 
discussed here for quantum dots is also closely related to 
the concept of doorway states in nuclear physics • 

We consider a dot-lead system described by the Hamil- 
tonian 



H = H, 



dot 
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Here, the dot is represented by the Hamiltonian 



A^dot 

Hdot = ^ (ei 

i=l 



fJ-)alai + C/ ^ alaia''-aj (2) 



in terms of the creation operators al of an electron in 
the i-th single-particle eigenstate of the dot with energy 
Ei, charging energy U — e^/C and chemical potential /i. 
For disordered dots its eigenstates and eigenvalues are 
usually obtained from a random matrix ensemble. The 
Hamiltonian of the fc-th lead reads 

oc oo 

j=i i=i 

where c*^^ is the creation operator of an electron on the 
j-th site of the fc-th ID lead, and t is the hopping matrix 
element in the lead. The coupling between the dot and 
the lead is contained in 



HhoupHng = E ^^.^^^t + h.C. 



(4) 



where the dot is assumed to be attached to the edge of the 
lead, and the coupling amplitude between the i-th orbital 
in the dot and the k-th lead is given by Vk^i- The N one- 
dimensional leads may also be connected by transverse 
hopping — i X)^/ c^^c*^"''^ + h.c. in order to turn 
them into a quasi-one dimensional lead with N channels. 

We begin with exact-diagonalization results for a non- 
interacting (U — 0) dot coupled to an external quasi- 
ID lead of varying width (i.e., number of channels N). 
Exact diagonalization can only treat finite systems, and 
therefore cannot deal with infinite leads. Nevertheless, as 
long as the level broadening in the dot is much larger than 
the level spacing in the lead, the description of the system 
is accurate. Diagonalizing of the Hamiltonian H gives 
the eigenvalues and eigenvectors |m) of the dot-lead 
system. The number of electrons on the dot at a given 
chemical potential is n = J^'^m^iJ^^Ji |(™|a|ai|m)p. 

In Fig. n we present the discrete increases An(/i, N) = 
J^^Q QQif (dn/dii') dii' in the occupation of a disordered 
quantum dot with A^dot = 16 orbitals (generated by a 
random matrix with a Gaussian distribution of width 
O.lt) as function of the chemical potential /i and N, the 
number of channels connected to the dot. In the absence 
of dot-lead coupling [N = 0), all increases An(^, A^) > 
occur at /i = e„i. The integer points on the x-axis corre- 
spond to N channels connected to the dot by couplings 
Vk,i drawn from a random Gaussian distribution with a 
zero mean and variance t. The non-integer values of N 
correspond to increasing logarithmically the values of the 
coupling of the iV -f 1 channel to the dot and its trans- 
verse hopping to the neighboring channel up to the full 
strength Vi_Ni+i,i and t at integer N. For integer values 
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FIG. 1: The increases An(/i, A*') in the occupation of a quan- 
tum dot with Adot ~ 16 orbitals as function of the chemical 
potential fi and the width A' of the quasi-one-dimensional 
lead (length M — 100) connected to the dot. Non-integer 
numbers of A^ correspond to having the connection of the last 
channel to the dot and its neighboring channel logarithmically 
increased (see text). The increase in the occupation is color 
coded as indicated in the figure. The symbols correspond to 
the derivative in the limit of V oo 



of A^ it can be clearly seen that there are Adot — A^ jumps 
in the occupation of the dot corresponding to the same 
number of states bound to the dot. As one couples an 
additional channel to the dot the energies of these states 
gradually change, until for some intermediate strength 
of coupling some of the states move abruptly, split, or 
disappear. At stronger couplings (i.e., close to the next 
integer A^) the ordered structure of the states reappears, 
with one less state than for the previous value of A'^. Sim- 
ilar behavior is seen when there is no transverse hopping 
between the channels (i.e., A^ independent ID leads). If 
Vk,i = t is independent of i or fc (i.e., all the couplings to 
the same lead or orbital are identical) a loss of a bound 
state occurs only when connecting the first lead to the 
dot. Attaching additional leads does not change the num- 
ber of bound states on the dot. 

To understand these numerical results we approach the 
system within the scattering-matrix formalism. For A^ 
propagating channels in the leads, the dot-lead system 
can be characterized by the N x N scattering matrix [l^ 
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dot 



(5) 



where V is an N x A^dot matrix describing the coupling 
of the channels to the Adot orbitals of the dot. The local 
DOS of the dot is given by dn/d^ ~ (l/7r)Imtr[^— i/dot + 
iTri'V^V]~^ . If, to be specific, we take i?dot to be diag- 
onal with random-matrix spectrum, the couplings Vk^i 
are essentially independent random Gaussian variables 
of variance v'^. We emphasize, however, that our results 
are not specific to a random-matrix spectrum. 
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For a weakly-coupled dot, we can compute the S- 
matrix by first diagonalizing the dot Hamiltonian i?dot- 
The broadening of the levels can then be obtained 
by treating tTTuV^V in first-order perturbation theory. 
To understand the opposite limit of strong coupling, 
we first diagonalize invV^V and subsequently account 
for ifdot perturbatively. We start by writing V'^ — 
(vi, V2, . . . , vjv) where the are A'dot-dimensional vec- 
tors. In obvious notation, we can then write iivV'^V — 
T^vY^f^-^ which shows that the A^dot x iVdot ma- 

trix TTvV^V has at most rank TV and generically only N 
nonzero eigenvalues ^ TrvN^otv'^ with i — 1, . . . ,N . 

Including i?dot perturbatively, we first need to diag- 
onalize i/dot in the (TVdot ~ A'^)-dimensional degenerate 
subspace of zero eigenvalues. In the limit A^dot ^ 
this leads to a random-matrix spectrum of A^dot — N real 
eigenvalues whose width and level spacing A equal 
those of the Hamiltonian i?dot of the uncoupled dot. 
Thus, in first-order perturbation theory we find TVdot — N 
infinitely sharp resonances in addition to N imaginary 
eigenvalues which lead to an extremely broad background 
(since A ^ Adot) in the local DOS. 

In second-order perturbation theory, these resonances 
acquire a width since iJdot couples the sharp resonances 
to the broad background. The resulting width can be 
easily estimated to be 



Ae 



N 

E 

1=1 



|(»|gdot|^)P 

-iXi 



N^A" 



(6) 



Here, we used that i/dot is a random-matrix Hamiltonian 
and defined the golden-rule width F = ttvNv^ of the 
eigenstates of the uncoupled dot. We assume strong 
coupling so that -kvv^ ^ A which allows us to neglect the 
unperturbed energy in the denominator. When this 
width remains small compared to the level spacing, i.e., 
when N^/ngdot ^ 1, we find N^ot—N isolated resonances 
in both conductance and local DOS even though dot and 
lead are very strongly coupled, in agreement with our 
numerical results. 

By an analogous argument one finds only one imagi- 
nary eigenvalue and hence A^dot — 1 sharp resonances in 
the non-generic cases in which the coupling Vk^i is inde- 
pendent of either channel k or dot level i. These res- 
onances have a width of order A/irgdof This width is 
smaller by a factor A^^ compared to the resonance width 
in the case of arbitrary dot-lead coupling. 

We now turn to the influence of the charging energy U. 
We compute the ground state for an interacting dot at- 
tached to several ID leads using an extension of a DMRG 
method developed for a single ID lead connected to a dot 
[l3| which will be described in detail elsewhere. As in the 
single lead case, the essence of the method is similar to 
the regular DMRG for ID systems Jjj. The main dif- 
ference is that in every iteration a site is added to each 
of the leads. Fig. [3 shows the occupation number n as 
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FIG. 2: The occupation n as function of the chemical po- 
tential for A^dot = 3 and N — 2. The dot levels have ener- 
gies £1 = —1.3t,e2 = —1.2t, and ea = —l.lt. The interaction 
U = 0.2t. In (a) the couplings are ,Vii = O.OSpt, V12 = 0.09pt, 
Vn = O.Olpt, V21 = O.OQpt, V22 = O.Olpt, V23 = O.Olpi, where 
p — 1, 2, 3 ... 12 for the different curves, while in (b) all the 
couplings are V = 0.05pt. 



function of fi for a A'dot = 3 dot attached to two leads 
for different values of the dot-lead coupling. 

In Fig. the case of non-identical couplings is pre- 
sented. The general behavior seen for the non-interacting 
case is repeated in the interacting case. For weak cou- 
pling, there are three discrete jumps in the occupation of 
the dot separated hy A + U. As the coupling increases 
only one discrete jump remains. Thus, the interactions 
in the dot do not eliminate the bound state. Moreover, 
interactions shift the position of the remaining jump to 
higher /i relative to the noninteracting case. 

If the couplings are symmetric to all leads the number 
of bound states is A^dot — 1 no matter how many leads 
are attached. This is illustrated in Fig. [2l3, where two 
discrete jumps remain even for strong coupling. These 
bound states are separated by a distance oi A + U as one 
expects from two bound states on an interacting dot. 
Similar Coulomb blockade behavior at strong coupling 
has been recently seen for the Kondo system [l5j . 

This behavior can be explained by extending the scat- 
tering theory above to include U within the Hartree ap- 
proximation. This approximation accounts for U by re- 
placing ei ^ ei + U'Yj_^^{bjbj) where bj annihilates an 
electron in the dot state Cj. The charging of the broad 
resonances can be neglected as long as we consider chem- 
ical potential changes which are small compared to the 
bandwidth. With this approximation, subsequent reso- 
nances are separated by A + U. 

Remarkably, the behavior changes for very large charg- 
ing energy U ^ Trz/Adot^^ where the superradiant state 
becomes Coulomb blockaded. In this regime, we observe 
numerically that the number of Coulomb-blockade steps 
equals the number of dot states A^dot, cf. Fig. 3(a). The 
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FIG. 3: Occupations vs. /i for A'dot = 3 (ei = — 1.41i,e2 = 
— lAt, 63 = — 1.39t) coupled to one lead (all couplings V — 
0.2t) with strong charging energy U = 0.6t. (a) Dot occupa- 
tion n = ni + n2 + ns. (b) Individual dot-orbital occupations 
ni (full line) , n2 (dashed line) , (long dashed) . (c) Occupa- 
tion in the strong-coupling basis: superradiant (full line) and 
subradiant states (dashed and long-dashed lines). 



width of the steps is of order Trz^iVdoti'^, large compared 
to the step widths of subradiant states at weaker charg- 
ing energy. Fig. 3(b) shows that the charging steps are 
due to simultaneous charging of all three dot orbitals 
which leads to an oscillatory structure in the occupations 
of the dot states. The origin of these oscillations can 
be traced by considering the occupations of the strong- 
coupling super- and subradiant states (i.e., the eigen- 
states of i'KvV'^V) as shown in Fig. 3(c). Clearly, the 
origin lies in oscillations in the occupation of the super- 
radiant state which, in this case, is a symmetric super- 
position of all dot orbitals. At the conductance step, one 
predominantly charges the superradiant state while in be- 
tween steps, there is a tendency to exchange occupations 
between the superradiant and a subradiant state. 

To understand this behavior, consider the configura- 
tion of the dot when /.t takes a value on the charging 
plateau n = \. In this case, the dot could either charge 
the superradiant or a subradiant state. In perturbation 
theory, the dot-lead coupling changes the (many-body) 
energy of these configurations (iSsupor and i?sub, respec- 
tively) due to virtually exciting an electron from the lead 
to the superradiant state, if the latter is unoccupied, or 
by virtually exciting an electron from the dot to the leads, 
if the superradiant state is occupied. (We neglect virtual 
processes involving subradiant states since their coupling 
to the leads is much weaker.) This gives an energy differ- 
ence = E,^^,, - E.^y, = f [ln(t/|e + ;7|) - ln(t/|e|)] 
between the two configurations, where e is the energy of 
the dot state relative to fj, and W is the width of the 
superradiant state. (The single-particle level spacing is 
neglected here.) Thus, for |e| < U/2 (|e| > U/2) occu- 
pying the superradiant (subradiant) state gives the lower 
energy so that the charging steps are due to charging of 



the superradiant state, while somewhere on the plateau 
the occupations of the superradiant and subradiant states 
are exchanged. The precise location of this switch is af- 
fected by the single-particle level spacing of the dot. This 
mechanism was considered by Silvestrov and Imry in a 
different context [l^ . 

We close with a discussion of possible experimental re- 
alizations. Generically, for a semiconductor quantum dot 
perfectly coupled to leads by quantum point contacts in 
the sense that there are N perfectly transmitting chan- 
nels, one finds that F ~ NA, leading to N^A/ttT ^ N. 
This shows that this situation is not in the strong cou- 
pling limit and thus, there are no sharp resonances in 
agreement with the description of this regime e.g. in Refs. 

One situation in which the effect discussed here can 
be observed is when, by a mesoscopic fluctuation or by 
symmetry, several dot levels bunch together so that their 
effective level spacing is much smaller than the average 
A. Alternatively, a situation with an anomalously small 
A can be engineered by a judicious choice of the device. 
E.g., one can think of a set of n quantum dots with weak 
interdot tunneling whose energy levels can be manipu- 
lated into almost degeneracy by a set of external gates. 
This realizes a situation with iVdot = n and strong cou- 
pling F ^ iV^A to the leads. Finally, one may also think 
of cases in which leads and "dot" are made from differ- 
ent materials, allowing for an independent manipulation 
of F and A. For example, when tunneling through a se- 
ries of identical impurities or a suitable molecule between 
metallic electrodes, one expects F to be enhanced by the 
large DOS in the metallic leads. 
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